12 Henon at tract or 



The chaotic phenomena of the Lorenz equations may be exhibited by even 
simpler systems. 

We now consider a discrete-time, 2-D mapping of the plane into itself. The 
points in M? are considered to be the the Poincare section of a flow in higher 
dimensions, say, R^. 

The restriction that > 2 for a strange attractor does not apply, because 
maps generate discrete points; thus the flow is not restricted by continuity 
(i.e., lines of points need not be parallel). 

12.1 The Henon map 

The discrete time, 2-D mapping of Henon is 

Xk+i = Yk + 1-aXl 

Yk+l = P^k 

• a controls the nonlinearity. 

• (3 controls the dissipation. 



Pictorially, we may consider a set of initial conditions given by an ellipse: 

Y 
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Now bend the elllipse, but preserve the area inside it (we shah soon quantify 
area preservation) : 



7' 

A 

Map Ti : X' = X 

Y' = l-aX'^ + Y 

Next, contract in the a;-direction < 1) 

Map T2 : X" = PX' 
Y" = Y' 

Finally, reorient along the X axis (i.e. flip across the diagonal). 

A 

Map T3 : X'" = Y" 
Y'" = X" 

The composite of these maps is 

T = T3oT2oTi. 

We readily find that T is the Henon map: 

X'" = l-aX^ + y 
Y'" = (3X 
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12.2 Dissipation 



The rate of dissipation may be quantified directly from the mapping via the 
Jacobian. 



We write the map as 



Xk+i = f{Xk,Yk) 
Yfc+i = 9{Xk,Yk) 



Infinitesimal changes in mapped quantities as a function of infinitesimal 
changes in inputs follow 

dXk dVk 

We may approximate, to first order, the increment AXk+i due to small in- 
crements (AXfc, AYfc) as 

df df 



dX, 



When {AXk, AY^) are perturbations about a point (xq, yo), we have, to first 
order. 



AXfc+i 




. An+i _ 





g'xS^o.yo) 5'k (^0,2/0) 



AXk 
AYk 



Rewrite simply as 



'Ax'' 




a b 




' Ax ' 


Ay' _ 




c d 




Ay _ 



Geometrically, this system describes the transformation of a rectangular area 
determined by the vertex (Ax, Ay) to a parallelogram as follows: 



Ay 



(Ax; Ay') 



(Ax Ay) 



Ax 
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Here we have taken account of transformations like 



(Ax,0) 
(0,A|/) 



{aAx, cAx 
(bAy, dAy) 



If the original rectangle has unit area (i.e., Ax Ay = 1), then the area of the 
parallelogram is given by the magnitude of the cross product of (a, c) and 
(6, d), or, in general, the Jacobian determinant 



Therefore 



J = 



\J\ > 1 
|J|<1 









dXk+i 


a b 




dXk 


dVk 


c d 






dVk+i 






dXk 





{xo,yo) 



^ dilation 



^ contraction 



For the Henon map, 

-2aXk 1 
P 

Thus areas are multipled at each iteration by 



J = 



After k iterations of the map, an initial area ao becomes 

ak = aQ\(3\^. 

12.3 Numerical simulations 



Henon chose a = 1.4, [3 = 0.3. The dissipation is thus considerably less than 
the factor of 10~^ in the Lorenz model. 



An illustration of the attractor is given by 



BPV, Figure VI. 19 
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Numerical simulations show the basin of attraction to be quite complex. 
Sensitivity to initial conditions is confirmed by 
BPV, Figure VI.20 

The weak dissipation allows one to see the fractal structure induced by the 
repetitive folding 

BPV, Figure VI.21 

Note the apparent scale-invariance: at each magnification of scale, we see 
that the upper line is composed of 3 separate lines. 

The fractal dimension D = 1.26. (We shall soon discuss how this is com- 
puted.) 

The action of the Henon map near the attractor is evident in the deformation 
of a small circle of initial conditions: 

BPV, Figure VI.22 

At the scale of the attractor we can see the combined effects of stretching and 
folding: 

BPV, Figure VI.23 
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